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ABSTRACT 

Context. Stardust is a class of presolar grains each of which presents an ideally uncontaminated stellar sample. Mainstream silicon 
carbide (SiC) Stardust formed in the extended envelopes of carbon-rich asymptotic giant branch (AGB) stars and incorporated the 
radioactive nucleus ^*A1 as a trace element. 

Aims. The aim of this paper is to analyse in detail the effect of nuclear uncertainties, in particular the large uncertainties of up to four 
orders of magnitude related to the ^^Alg(p, 7)^'Si reaction rate, on the production of Al in AGB stars and compare model predictions 
to data obtained from laboratory analysis of SiC Stardust grains. Stellar uncertainties are also briefly discussed. 
Methods. We use a detailed nucleosynthesis postprocessing code to calculate the ^'A1/-^A1 ratios at the surface of AGB stars of 
different masses (M = 1.75, 3, and 5 Mq) and metallicities (Z = 0.02, 0.012, and 0.008). 

Results. For the lower limit and recommended value of the ^^Alg(p,y)-^Si reaction rate, the predicted ^*A1/^'A1 ratios replicate the 
upper values of the range of the A1/-'A1 ratios measured in SiC grains. For the upper limit of the ^^Alg(p, 7)^'Si reaction rate, instead, 
the predicted '*A1/^'A1 ratios are x 100 times lower and lie below the range observed in SiC grains. When considering models of 
different masses and metallicities, the spread of more than an order of magnitude in the -*A1/^' Al ratios measured in stellar SiC grains 
is not reproduced. 

Conclusions. We propose two scenarios to explain the spread of the ^'A1/-^A1 ratios observed in mainstream SiC, depending on the 
choice of the ^^Alg+p reaction rate. One involves different times of Stardust formation, the other involves extra-mixing processes. 
Stronger conclusions on the interpretation of the Al composition of AGB Stardust will be possible after more information is available 
from future nuclear experiments on the ^^Alg+p reaction. 

Key words, nuclear reactions, nucleosynthesis, abundances - stars: AGB and post-AGB 



1. Introduction 

Meteoritic stellar grains are solid samples of stars that can be 
studied in terrestrial laboratories. Condensed in the cooling gas 
outflows from ancient stars, they became part of the interstellar 
medium from which the Solar System formed some 4.6 billion 
years ago. Because they were encapsuled in primitive meteorites 
they remained ideally uncontaminated by Solar System material. 
The highly unusual isotopic ratios, with respect to solar, found 
in these "Stardust" grains indicate that they are of stellar ori- 
gin and can therefore be used as a diagnostic tool for verifying 
predictions from models of stellar evolution and nucleosynthe- 
sis (c f. [Anders & ZinneBEMU IZinneij|l998t IClavton & NittieJ 
l2004tlLugaroii2005 ). Stardust grains come in many mineralog- 
ical flavours but of particular interest in this work are silicon 
carbide (SiC) grains. 

"Mainstream" SiC grains (>90% of stellar SiC grains) con- 
tain isotopic abundances of heavy elements characteristic of the 
slow neutron capture process (the s process) and are thus be- 
lieved tojiavejjrminated from as ymptotic giant branch (AGB) 
stars dLugaro et alj 1 19991 l2003ah . which show enrichments at 
their surface of jf-process elements such as Zr, Ba, and even the 
unstable Tc dMerilll 119521: ISmifli & Lamber3ll990t iBusso et alj 



I2001h . These stars are evol ved giants in t he final nuclear burn- 
ing stage of evolution (see lHerwidl2005L for a recent review). 
Briefly, the AGB phase is important because of the occurrence 
of instabilities of the He-burning shell, known as thermal pulses 
(TPs). Inbetween TPs the H-shell provides most of the stellar 
luminosity. After the occurrence of a thermal pulse the third 
dredge-up (TDU) may occur, where the products of the partial 
He-burning such as '^C, along with other nucleosynthetic prod- 
ucts such as i-process elements, are mixed to the stellar surface. 

An interesting nucleosynthetic product of AGB stars is the 
radioactive isotope ^^Al (T1/2 = 0.717 Myr). ^''Al is produced 
by proton captures on ^^Mg during H burning, when it can 
also be consumed by proton captures, depending on the rate 
of the ^^Al+p reaction, and destroyed by neutron captures in 
the thermal pulse because of the relatively high ^^Al(n,p)^^Mg 
and ^^Al(n,Q')^^Na rates. Neutrons in the TP are provided by the 
^^Ne(Q',n)^^Mg reaction if the temperature exceeds « 300 mil- 
lion K. When the TDU occurs the Al in the thin top layer of the 
intershell region not involved in the convective pulse, together 
with the ^^Al that survived neutron captures in the TP, is carried 
to the surface. For a detailed description of the nucleosynthesis 
of^^Al in AGB stars see Mpwlavi & Meynet (2000). The combi- 
nation of nucleosynthesis in the intershell and the occuiTence of 
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Murchison : 71 
Orgueil : 57 
Both : 128 




Fig. 1, Histogram of binned ^*A1/^^ Al ratios measured in main- 
stream SiC Stardust grains found in the Murchison and Orgueil 
meteorites. The total number of grains for which the ^''Al/^^Al 
ratios were measured is indicated in the upper right hand corner. 



the TDU allows AGB stars of masses roughly between 1.5 M^T) 
and 4 Mq to eventually become carbon rich (Groenewege n et al] 
Il995h . which is a necessary condition for the formation of SiC 
grains. Al is incorporated in SiC grains as a trace element. 

Since the abundance of Mg in SiC grains is much lower than 
that of Al, excesses in ^*Mg, together with solar ^^Mg/^"^Mg ra- 
tios, observed in meteoritic Stardust SiC grains are a measure 
of the ^^Al abundance at the time and place of formation of the 
grain. Note that if spallation reactions on the grains during their 
residence time in the interstellar medium had been responsible 
for the observed huge excesses in ^^Mg, this would also result in 
variations of the ^^Mg/^'*Mg ratios, which are not observed. The 
distrib ution of the available ^^Al/^^Al data for mainstream SiC 
grains dHoppe et alJll994l : lHuss et alJI 19971 ) is presented in Fig. 
[1] The majority of the grains are distributed between ratios of x 
10""* to 2 X 10"^. The ^^Al/^^Al ratio has also recently been mea- 
sured in a few SiC grains of the rare ty pe Z, which are bel ieved 
to have originated in low-Z AGB stars (IZinner et al.ll2007l) . The 
data points cover the same spread as the mainstream grains and 
low metallicity AGB models give similar results to our x solar 
metallicity models. Hence, our discussion can be applied to Z 
grains as well. 

The aim of this paper is to analyse in detail the effect of 
the uncertainty in the ^^Mg(p,y)^^Al and ^^Alg(p,y)^^Si reac- 
tion rates on the production of ^''Al in carbon rich AGB stars 
and compare the model predictions to the data obtained from 
mainstream SiC grains. While the -^^Mg+p is uncertain by a fac- 
tor of Si 2 at the temperature of interest for H burning during 
the AGB stage 60 million K), the ^^Al+p reaction rate is un- 
certain by four orders of magnitude (see Fig. |2]i due to possi- 
ble contributions from as yet unobserve d low-energy resonances 
dlliadis et alj|2001l: lAngulo et al.lll999l) . We evaluate the effect 
of such huge uncertainties on the outputs of AGB models and 
see if any constraints can be derived from comparison of these 
models to the SiC data. Stellar model uncertainties will also be 
briefly discussed. 

The study of radioactive isotopes in Stardust grains is 
of interest also because they represent clocks to measure 
the timescale for the formation of dust around AGB stars 
dZinner et al ] |2006al:lDavis & Gamnoir2006 l). It is still very diffi- 
cult to pin down the exact mechanism by which dust is formed 



around AGB stars, see for example discussion in iNuth et all 
(|200l),and any insight on the timescale of grain formation there- 
fore represents a useful constraint. 

We also note that ^*A1 is produced by proton captures at the 
base of the envelope (hot bottom burning, HBB) in intermedi- 
ate mass AGB stars (M ^ 5 Mq) and in AGB stars of lower 
masses if extra mixing is invoked, which is usually done to ex- 
plain the composition of a particular fraction of Stardust ox- 
ide g rains showing a strong depletion of ^^ O along with high 
^^Al dNittler et al.1 119971: iNollett et al.ll2003l) . Here we present 
only one model of an intermediate mass AGB star because they 
are not main producers of SiC Stardust. In fact, HBB prevents 
massive, metal-rich (Z > 0.004) AGB stars from developing a 
carbon rich envelope by converting C into N. Moreover, the 
isotopic signatures of C, N, Si, and heavy elements in main- 
stream Si C grains cannot be reconcil ed with massive AGB par- 
ent stars dLugaro et al.lll999l l2003al) . For example, HBB pro- 
duces '^C/'^C in the range 3 to 10, while the mainstream SiC 
range is 20 to 100 and high neutron densities in the thermal 
pulses produce enhancements in the neutron-rich isotope "^^Zr, 
which is instead observed to be depleted in mainstream SiC 
grains. We have not included extra mixing in our models, but 
will d iscuss its pos sible implications in §4, based on the mod- 
els of 'N olle tt et all d20b3). In §2 we present the details of the 
methods and models we have used. In §3 we show the results we 
obtained, which we discuss in §4. We conclude with §5 where 
we present a summary and our main conclusions. 



2. Methods and Models 

2.1. Evolutionary and Nucleosynthesis codes 

We calculate the nucleosynthesis with a detailed post-processing 
code for which the stellar structure inputs were calculated be- 
forehand. The stellar structure program we use has a small net- 
work of 6 species: H, ^He, "^He, C, N, and O (see e.g. iLattanziol 
il986i) . involved in the main energy generation. To compute 
abundances of more species we use a post-processi ng code. The 
code is described in detail in lLattanzio et al.l (Il996l) . Briefly: we 
input the structure (e.g. temperature, density, details of the con- 
vective regions, position of the H and He-burning shells as a 
function of interior mass and time) from the stellar evolution 
code to compute the abundances of species not involved in en- 
ergy generating reactions. The post-processing step computes its 
own mass mesh with sufficient resolution in each burning shell 
(around 25 mass shells) to adequately resolve changes in abun- 
dances owing to nuclear burning. Because we assume the struc- 
ture is fixed in the post-processing step we assume that the ex- 
tra species and reaction rates added do not change the structure. 
For this reason we concentrate on uncertainties for reactions that 
produce negligible energy, such as those involved in the NeNa 
and MgAl chains, and neutron capture rea ctions. The details o f 
the nucleosynthesis network are outlined in lLugaro et al.l d2004l) . 
however, we remind the reader that we include 59 light nuclei 
and 14 iron-group species. We also add the fictional particle g 
to count the number of neutron captures occurring beyond ^^Ni 
(Lat tanzio et a l. 1996; Luga ro et al.ll2004l) in a similar manner to 
Jorissen & Arnould . C1989) . 

The bulk of our 527 r e action rates are from the REACLIB 
tables (Thielemann et al.' 1986) updated as described in 
Lugaro et al. (2004) and Karakas et al. (2006). 

In our calculations we used models of 1.75 Mq with a 
metaflicity (Z) of 0.008, of 3 Mq with Z = 0.02, 0.012, and 
0.008 and of 5 Mq with Z = 0.02 and the mass loss pre- 
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scription from I Vassiliadis & WoodI ( Il993h . The 1.75 Mq and 3 
Mq models produce a carbon over oxygen ratio in excess of 
unity. More informat i on on these stel l ar mo dels c an be found in 
iLugaro et al] (l2003bh . iKarakas et all (l2006 l). and ' Karakas et all 
(l2007l) . Models of low-mass AGB stars of approximately solar 
metallicity are shown to be the best to reproduce various features 
of Star dust mainstream SiC grains: from the He and Ne compo- 
sition dGalUno etal.lll990h . to the '^C/'^C ratios, which are the 
same as observed in Gala ctic C stars with metallicity close to 
solar (Hoppe & Ot3ll997l) . to the heavy element compositions 
dLu^aro et al.l2003 a'). Consequently, we focus on our AGB mod- 
els of low mass and Z close to solar that beco me carbon rich. The 
initial abundances we use were taken from lAnders & Grevess^ 
(119891) for Z = 0.02 and we assume scaled solar fo r the Z ^ 0.008 
model s. Th e Z - 0.0 1 2 mo dels are based on the lAsplund et al.l 
(l2005h and iLoddersI (l2003l) solar abundance table which pre- 
scribes a solar metallicity of Z = 0.012. 

2.2. The reaction rates 



The 



rates of interest are those of the ^^Mg(p,'y)^^Alg, 
^^Mg(/?,y)^''Al„„ and ^^Al^(/7, y)2^Al reactions. ^''Al^, is the 
ground state of ^^Al, whereas ^*A1,„ is the metastable state with 
a half life of 6.3452 s. Since the metastable state of ^^Al is very 
short lived, the ^^Mg{p, y)^^Alm rate results in the production of 
essentially no ^''Al, but rather ^^Mg. Whenever we use ^^Al with 
no subscript we mean the total sum of the metastable and ground 



state of Al. In practice this means Alj, because the isomeric 
state i s very unstable. The rates were taken from ilhadis et"aLl 
(120011) . We calculated models using all the nine combinations 
of the lower limits (LL), recommended values (RC), and up- 
per limits (UL) of the rates, where the rate errors of the reac- 
tions ^-''Mg(p,7)^^Al^ and ^-''Mg(p,7)^^Al„, are correlated. This 
assumption is justified since the uncertainties arise from the en- 
trance channel partial width which is common to both interac- 
tions. The bulk of the rates in our code are calculated using fits in 
the 7-coefficient format of REACLIB; the rates of interest were 
read directly from the tabulated rates. 

Figure |2] shows the LL, RC, and UL rates for the 
^^Mgip, y)^^Alg and ^^Alg(p, 7')^^Si reactions for the tempera- 
ture range dominant in the H burning shell during the AGB 
phase, where 6x10^ K is approximately the peak temperature 
for all of the used models, except the one with 5 Mq and Z - 
0.02 where it is 8x10^ K. The upper Hmit of the ^^Alg{p, y)2^Si 
rate is up to four orders of magnitude higher than the lower 
limit and results, as shown below, in a two orders of magnitude 
lower abundance of ^^Al. These large uncertainties arise from as 
yet undetected low-energy resonances. In particular, an expected 
resonance at a center-of-mass energy of 94 ke V, with a predicted 
upper limit of its strength of a>y <10"'* eV, seems to play the 
most important role. 
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Fig. 2. The lower limit (bottom line), recommended value (mid- 
dle line), and upper limit (upper line) of the ^^Mg(/?, y)^^Alg 
(dashed lines) and ^*Alg(/?,7)^^Si (solid lines) reaction rates. 
The vertical line denotes the approximate temperature of interest 
(T=*6xlO^K). 




C/0 

Fig. 3. Surface ^^Al/^^Al ratio versus C/O ratios for a star of 
mass 3 Mq, metallicity Z - 0.02, and a partial mixing zone 
(PMZ) of IxlO"-' (see 33.1b using all combinations of upper, 
lower and recommended values for the rates under considera- 
tion. At the end of some selected lines a label indicates the val- 
ues used for the ^^Mg(/?, y)^^Als, and Ali,(/?, y)^^Si rates, in that 
order, where LL = lower limit, RC - recommended, and UL = 
upper limit. The grey crosshatched box is a schematic represen- 
tation of the range of the ^^Al/^^Al values observed in Stardust 
maintream SiC grains (see Fig. [T]l, which can form only when 
C/0>1. 



3. Results 

Figure [3] shows our predictions for the ^^Al/^^Al ratio at the sur- 
face of the star as a function of C/O for the 3 Mq Z-0.02 model 
and all nine different combinations of the LL, RC, and UL of the 
rates. All the lines show a very similar trend: the ^^Al/^^Al ra- 
tio increases sharply with the first few TDUs and then becomes 
approximately constant. This is because during the first few TPs 
the mass fraction of ^*A1 dredged up to the surface is 6.2x10"^ 
in the small region (^ 10 Mq) of the H-burning ashes (region 
A of iMowlavi & Mevnetll2000) and ^ 1.8x10"^ in the rest of 



the intersh ell (region D of iMowlavi & Mevnelll2000l) . (See also 
Table 2 of iLugaroet al.l200l l). However, as the pulse number in- 
creases, the temperature at the base of the convective intershell 
region increases, and the ^^Ne(Q',n)^^Mg reaction becomes more 
active, this frees up more neutrons resulting in more ^^Al de- 
struction by neutron capture. The abundance of ^^Al in the last 
computed TP is in fact ^ 2x10 in the intershell. The TDU of 
this small amount of ^^Al is just enough to ensure that the small 
fraction (» 5%) of ^^Al at the stellar surface that decays during 
the interpulse period is replenished. As a result the prediction 
lines flatten. 
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Fig. 4. ^^Al/^^ Al abundance ratios for models of different mass 
and metallicity. Two calculations are shown for each model: the 
top line represents the reaction rate combination: UL LL and the 
bottom line the reaction rate combination: LL UL (see Fig. [3]). 
The fact that the 3 Mq Z-Q.Q2 line is slighlty higher than shown 
in Fig.[3]is due to the fact that a partial mixing zone was included 
in the calculations shown in Fig. [3] while no partial mixing zone 
is included in the calculations shown in this figure. 



Two main classes of prediction lines can be distinguished: 

1. When considering the models computed using the RC and 
LL values for the ^^Alg(p, 7')^^Si rate the variation among 
these is of a factor of a mere 1.32. Of this, a factor of ~ 1.26 
is derived when varying the ^^Mg(p, 7)^^Ali, reaction rate be- 
tween the LL and the UL, while a factor of ^ 1.05 is derived 
when varying the ^^Alg{p, y)^^Si reaction rate between the 
LL and the RC rate. All the lines lie around the upper end of 
the SiC grain data range. 

2. When considering the models computed using the UL values 
for the ^^Algip, 'y)^^Si rate, the ^^Al/^^Al ratio is roughly two 
orders of magnitude smaller than the ratios computed using 
the RC and the LL values of the rate, hence falling below the 
range of values observed in mainstream SiC grains. 

In Fig. |4] we present the results for different masses and 
metallicities. We can see the same behaviour here as described 
above: the upper limit of the ^^Alg{p, y)^'Si reaction rate always 
results in a factor of w 100 less ^^Al, except for the 5 Mq model, 
in which case the difference is of a factor of « 10. Note that the 
production of -^Al in this massive AGB model is due to the oper- 
ation of the second and third dredge-up and not due to hot bottom 
burning because the temperature at the bottom of the convective 
envelope is not hig h enough to produce ^^Al (see discussion in 
iLugaro et al.l2007l) . This explains why our detailed model shows 
a large effect in the ^^Al abundance as a result of the uncertain- 
ties in the^^Alg(p, jj^^S i reaction rate, while the synthetic mod- 
els of llzzard et al.l ( 2007b . which only describe the effect of hot 
bottom burning, do not show this effect. 

All the models computed using the RC or LL for the 
^^Algip, 7)^^81 reaction rate produce ^^Al/^^Al ratios at the up- 
per end of the observed distri bution. Similar resu lt s wer e ob- 
tained (for t he RC case) by iMowlavi & Mevnetl ( |2000|) and 
IZinner et al.l (!200 7|) for s everal models of different masses and 
Z, and bv .Cristalloet al.l(l2006h for a 2 M© Z = 0.015 model. 



3. 1 . Other uncertainties 

Some uncertainty is introduced by the neutron capture re- 
action rates. We use for these rat e s the estimates of 
iKoehler et all d 19971) and ISkelton et al] (Il987h . These works 
suggest for th e ^^A le (», p) ^^Mg rate a value similar to that of 
[Caughlan & F owled ([T988I) . ^ 250 mbarn (at 23 keV, taken 
as typical temperature in the thermal pulse) and for the 
^^Al j,(«, Q-)^-^Na rate a value r oughly a factor of two higher 
than lCaughlan & Fowled ( Il988l) . ^ 180 mbarn (at 23 keV). The 
^*Aly(n, y)^^Al rate is much smaller than the aforementioned 
neutron capture channels: w 4.5 mbarn (at 23 keV, Bao et al] 
:2000). We conservatively estimated the uncertainties of these 
rates to be of a factor of two both above and below the values 
we use. When we multiplied the rate of all n-capture rates on 
^^Al by a factor of two we ended up with half as much ^^Al, 
similarly when we decreased the «-capture rates by a factor of 
two, we ended up with twice as much ^^Al. Overall, the uncer- 
tainty is at most of a factor of four. 

Since neutrons released by the ^^Ne(Q;,n)^^Mg reaction de- 
stroy ^^Al, the uncertainty of this reaction rate introduces an- 
other uncertainty in the ^^Al abundance. B etween the lower and 
upper limit (a range of about a factor of 15, iKarakas et alj|2006l) 
of this rate the Al abundance varies with a factor of 1 .8 for the 
3 Mq and Z = 0.02 model. 

Note that the uncertainties of a factor of ^ 20 (llliadis et al.l 

1200 Ih in the ^''Al+p reaction rate do not change the ^^Al abun- 
dance because this reaction is too slow in the H-burning shell of 
AGB stars. 

As for stellar uncertainties: we artificially included a partial 
mixing zone in th e top layers of the intershell in the way de- 
scribed in detail bv lLugaro et alj (12004}) . There, protons combine 
with '^C to make '■'C, which releases neutrons for the s process 
via '^C(Q',n)'*0 during the interpulse periods, however, the na- 
ture of the mixing is still debated (Herwig 2005). Varying the 
size of the '^C pocket from to 2x10"^ introduces a spread in 
the Al/^^ Al ratio of at most a factor of 1 .4. Varying the proton 
profile in the partial mixing zone changes the relative importance 
of the '"^N-p oor and '^N-rich re gions of the pocket. As shown 
in Fig. 1 of Gorielv & Mowlavil Q.000), in the '"^N-poor region 
^^Al is completely destroyed, while in the '"^N-rich region it is 
destroyed to an abundance of the order of 10"^ in number. Thus, 
changes in the proton profile would not have a significant impact 
on the overall ^^Al abundance in these stars. 

Mass loss and third dredge-up are very uncertain physical 
features of AGB stars. However, we expect that changes in the 
TDU efficiency, caused by variations in the input physics, and 
changes of the mass loss values will not affect our results sig- 
nificantly. This is because we are looking at C-rich stars and the 
C/0>1 constraint sets the dilution factor in our models to 1 
part of intershell material to ^ 30 parts of envelope material. This 
dilution factor would necessarily produce ^^Al/^^Al ^ 10"^, in- 
depently of which TDU efficiency and mass loss rate have been 
employed to achieve it. 

Extra-mixing processes due to rotation or other mechanisms 
may also be at work in these stars. Extra mixing is the hypothesis 
that some of the material from the convective envelope is mixed 
into the radiative layer that resides on top of the hydrogen shell. 
It is also commonly referred to as: "deep mixing" and "cool bot- 
tom processing". Extra mixing was originally introduced into the 
first giant phase of evolution to explain several abundance pecu- 
liarities, including lower than predicted '^C/'^C rati os in first gi- 
ant branch stars (see e.g. lSweigart & Mengel|[T979l) . Extra mix- 
ing has been invoked for AGB stars to explain the ^^C/'^C ratio 
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in C-rich stars ('Abia & Isern'2000'), as well as in mainstream and 
Z SiC grains (Zinner et al. 2006b), the O composition of AGB 
stars (see e.g lWasserburg et al.lll995h and t he O and Al compo - 
sition of a fraction of stellar oxide grains (iNoUettet al.ll2003h . 
The possible effect of extra-mixing on the ^^Al abundance will 
be discussed in §4.2. 



4. Discussion 

With the exception of the 5 Mq model, the models presented 
in Fig. m should represent the set of stars, in terms of mass 
and metallicity, responsible for mainstream SiC Stardust grains. 
However, the spread of the -^Al/^^Al SiC grain data is much 
larger than that of the models. The spread of the data is not an 
error range, but arises from considering the large number of sin- 
gle grains (128, to be exact, see Fig.[T]). The ^^Al/^^Al abundance 
ratios range a factor of ^ 20 in the data, whereas the models show 
at best a spread of about a factor of 2. We propose two scenarios 
to explain the observed range: 

4.1. The ^''Alg+p reaction rate corresponds to its current LL 
or RC value 

In this case we need to explain the fact that the data extends 
below the prediction lines, which may reflect the fact that star- 
dust SiC grains have a long interval of formation time between 
a time close to when TDU has enriched the stellar winds with 
^^Al and approximately two million years later This time delay 
would allow the decay of ^^Al into ~^Mg, the latter of which is 
not incorporated in the grains. Since the mechanism by which 
the large (up to 25 /vm) Stardust SiC grains found in meteorites 
are formed around single AGB stars is not understood, and in any 
case it probably does not involve timescales longer than ~ 10^ 
yr (fNuth et al. 2006), we speculate upon different possibilities. 

A long interval of formation time could be achieved if the 
grains were formed either (i) in the winds of extrinsic carbon 
stars formed in a binary system (in this case the absence of ^^Al 
would be conceptually equivalent to that of Tc in these stars) or, 
perhaps, (ii) in a long-lived circumbinary disk that formed af- 
ter one of the binary stars evolved through the AGB phase. Note 
that circumbinary disks have already been proposed as the site of 
origin of at least some meteoritic Stardust i n order to explain the 
large observed sizes of the grains > 0.5 yum ( lJurall997 ).The long 
timescale of Stardust grain formation proposed here is not ruled 
out by the analysis of radioactive heavy nucl ei in mainstream 
SiC, w hich have been discussed in detail by iDavis & Gallinol 
(l2006h . 

4.2. The ^^Aig(p,yf'^ Si reaction rate corresponds to its 
current UL 

In this case we need to explain the fact that the data extends 
above the prediction lines, which may reflect the occurrence of 
extra-mixing processes. Extra mixing could produce ^^Al/^^Al 
ratios higher than our predictions and the spread in the data could 
be achieved by a range of temperatures (i.e. depths) to which 
the mixed material is exposed. The main problem with invoking 
extra mixing is that we still miss the physical mechanism behind 
the process, even if some steps forwa rd in this searc h have been 
made recently (Eggleton et al. 2006; Charbonnel & Zahn 2007; 
iBusso et al.ll2007l)" ." 

Note that assuming the UL for the ^^Alg{p,j)^^Si reaction 
rate does not affect Al production in supernovae, since most of 
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Fig. 5. Upper limit, recommended value, and lower limit of the 
^*Alg(/?, y)^^Si reaction rate (dashed lines) as in Fig. ID with 10 
logarithmically equidistant steps in between the LL and the UL 
(solid lines). The shaded areas represent roughly which value of 
the ^^Alj(p, 'x)^^Si reaction rate corresponds to which scenario. 
The vertical yellow line denotes the approximate temperature of 
interest (as in Fig|2|l. 

this nucleus is produced during the explosion where the main de- 
struction channel for ^^Al is neutron capture dLimongi & Chieffil 
i2006) . Thus, the present scenario would not be in cont rast 
with the match to the observed live ^^Al in the Galaxy (Diehl| 
l2006h. On the other hand, production of ^''Al in Wolf-Rayet stars 
jAmould et al.>2006) and in intermediate-mass AGB stars by hot 
bottom b urning would be aff ected (for a detailed analysis on the 
latter see llzzard et"aLll2007h . As a consequence, the discussion 
on t he origin of Stardust spinel grain OC2 from a massive AGB 
star dLugaro et al.ll2007h would have to be revised. 

5. Summary and conclusions 

We have shown that the large uncertainty of the ^^Ali,-!-/? reaction 
rate has important implications when comparing AGB models to 
the ^*A1/^^A1 ratios observed in Stardust mainstream SiC grains. 
We have presented two scenarios to explain the observed distri- 
bution: one involves using the LL or the RC rate and invoking a 
relatively long timescale of grain formation, possibly connected 
to processes occurring in binary systems. The other involves us- 
ing the UL of the rate and invoking extra-mixing processes. 

To check in details which scenario would correspond to 
which values of the rate we have computed the ^''Al/^^Al ratios 
in the 3 Mq, Z =0.02 model for ten different rates calculated as 
logarithmically equidistant steps in between the LL and UL. The 
results are shown in Fig.|5] The decay scenario has to be invoked 
with rate values ranging from the LL up to the 6''' intermediate 
line from below. For rates from roughly the 8''' intermediate line 
up to the UL the ^''Al/^^Al ratio lies below the SiC data, and 
hence would imply the extra mixing scenario. 

Our work calls for more laboratory measurements of the 
^^Al/^^Al ratios in Stardust SiC grains of different sizes and for 
improved measurements of the ^^Alg+p reaction rates. A di- 
rect measurement of the 94 keV resonance may prove challeng- 
ing, considering the difficult preparation of a radioactive ^^Al 
beam-stop target that will not degrade under intense proton bom- 
bardment. An interesting attempt at studying the nuclear struc- 
ture of the corresponding compound level in ^^Si by using the 
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^^Alp( ^He,ij)^^Si reaction has been reported in IVogelaar et alj 
(Il996l) . Unfortunately, the contamination of their evaporated 
transmission target precluded the observation of transfers to any 
threshold states. A remeasurement of the transfer reaction by us- 
ing an isotopically and elementally pure implanted ^''Al target 
seems more promising. 

Once the estimate of the ^^Algip,-/) reaction is less uncertain, 
we will be able to readdress the issue of Al production in AGB 
stars and deliver stronger conclusions on the implications for 
dust formation around AGB stars and extra-mixing processes. 

Acknowledgements. ML gratefully acknowledges the support of NWO through 
the VENI grant. We thank Peter Hoppe and Gary Huss for providing us with the 
datatables of the grain measurements and for helpful and instructive discussion. 
We thank Alessandro Chieffi and Marco Limongi for discussion. AIK gratefully 
acknowledges the financial support of the Australian Research Council from the 
Discovery Project funding scheme (project number DP()664105). We thank the 
referee. Marcel Amould, for a detailed report which helped to improve the paper. 



References 

Abia, C, & Isern, J. 2000, lAUS. 177. 89A 

Anders, E., & Grevesse, N. 1989, Geochim. Cosmochim. Acta, 53, 197 

Anders, E., & Zinner, E. 1993, Meteoritics, 28, 490 

Angulo, C, Amould, M. Rayet, M. et al. 1999. Nucl. Phys. A, 656, 3 

Arnould, M., Goriely, S., & Meynet, G. 2006, A&A, 453, 653 

Asplund, M., Grevesse, N. & Sauval, A. J. 2005, in Cosmic Abundances as 

Records of Stellar Evolution and Nucleosynthesis in honor of David L. 

Lambert, eds. T. G. Barnes IH, & F. N. Bash (San Francisco: Astronomical 

Society of the Pacific) 336, 25 
Bao, Z. Y., Beer, H., Kappeler, F. et al. 2000, Atomic Data and Nuclear Data 

Tables, 76, 70 

Busso, M., Gallino, R., & Wasserburg, G. J. 1999, ARA&A, 37, 239 
Busso, M., Gallino, R., Lambert, D. L., Travaglio, C, & Smith, V. V. 2001, ApJ, 
557, 802 

Busso, M., Wasserburg, G. J., Nollett, K. M., & Calandra, A. 2007, ApJ, accepted 

astro-ph:0708.2949 
Cannon, R.C. 1993, MNRAS, 263, 817 

Caughlan, G.R., & Fowler, W.A. 1988, Atom. Data Nucl. Data Tables, 40, 283 
Charbonnel, C, & Zahn, J.-R 1007, A&A, in press |astro-ph/0703302 1 
Clayton, D. D., & Nitder, L. R. 2004, ARA&A, 42, 39 

Ciistallo, S., Galhno, R., Straniero, O., Piersanti, L., & Domuiguez, I. 2006, 

MemSAIt, 77, 774 
Davis, A. M., & Galhno, R. 2006, MemSAIt, 77, 885 
Diehl, R. 2006, New Astronomy Reviews, 50, 534 

Eggleton, P R, Dearborn, D. S. R, & Lattanzio, J.C. 2006, Science, 314, 1580 
Frost, C. A., & Lattanzio, J. C. 1996, ApJ, 473, 383 

Gallino, R., Busso, M., Picchio, G., & Raiteri, C. M. 1990, Nature 348, 298 
Groenewegen, M. A. T, van den Hoek, L. B., & de Jong, T. 1995, A&A, 293, 
381 

Goriely, S. & Mowlavi, N. 2000, A&A, 362, 599 
Herwig, F 2005, ARA&A, 43, 435 

Hoppe, R, Amari, S., Zinner, E., Ireland, T, & Lewis, R. S. 1994, ApJ, 430, 870 
Hoppe, P., & Ott, U. 1997, in "Astrophysical Implications of the Laboratory 

Study of Presolar Materials", ed. Thomas J. Bernatowicz & Ernst Zinner, 

American Institute of Physics Conference proceedings 402, 27 
Hoppe, P, Strebel, R., Eberhardt, P, Amari, S., & Lewis, R.S. 1996, Geochim. 

Cosmochim. Acta, 60, 883 
Huss, G.R., Hutcheon. I.D., & G. J. Wasserburg, G. J. 1997, Geochim. 

Cosmochim. Acta, 61, 5117 
lUadis, C, D'Auria, J. M., Starrfield, S., Thompson, W. J., & Wiescher, M. 2001, 

ApJS, 134, 151 

Izzard, R. G., Lugaro, M., Karakas, A. I., Ihadis, C, & van Raai, M. 2007, A&A, 
466, 641 

Jorissen A., & Amould, M. 1989. A&A, 221, 161 

Jura, M. 1997, in "Astrophysical Implications of the Laboratory Study of 
Presolar Materials", ed. Thomas J. Bernatowicz & Ernst Zinner, American 
Institute of Physics Conference proceedings 402, 379 
Karakas, A. I.. Lugaro, M., & Gallino, R. 2007, ApJL, 656, 73 
Karakas, A. I., Lugaro, M., Wiescher, M., Goeires, J., & Ugalde, C. 2006, ApJ, 
643,471 

Koehler, PE., Kavanagh, R.W., & Vogelaar, R.B. et al. 1997, Phys. Rev. C, 56, 
1138 

Lattanzio, J. 1986, ApJ, 311, 708 



Lattanzio, J., Frost, C, Cannon, R., & Wood, P. R. 1996, Memorie delta Societa 

Astronomia Italiana, 67, 729 
Limongi. M., & Chiefii, A. 2006, ApJ, 647, 483 
Lodders, K. 2003, ApJ, 591, 1220 

Lugaro, M. 2005, Stardust from Meteorites: An Introduction to Presolar Grains, 

Series in Astronomy and Astrophysics, Vol. 9 (World Scientific, Singapore) 
Lugaro, M., Karakas, A., Champagne, A., Lattanzio J., & Cannon, R. 2001, 

MemSAIt, 72, 319 
Lugaro, M., Karakas, A. I., Nitder, L. R. et al. 2007, A&A, 461, 657 
Lugaro, M., Davis, A. M., Gallino, R. et al. 2003a, ApJ, 593, 486 
Lugaro, M., Herwig, F, Lattanzio, J. C, Gallino, R., & Straniero, O. 2003b, ApJ, 
586, 1305 

Lugaro, M., Ugalde, C, Karakas, A. I. et al. 2004b, ApJ. 615, 934 
Lugaro, M., Zinner, E., Gallino, R., & Amari, S. 1999, ApJ, 527, 369 
Men-ill, R W. 1952, Science, 1 15, 484 
Mowlavi, N., & Meynet, G. 2000, A&A, 361, 959 

Nittler, L. R., Alexander, C. M.O'D., Gao, X.. Walker, R. M., & Zinner, E. 1997, 
ApJ, 483, 475 

Nollett, K. M., Busso, M., & Wasserburg, G. J. 2003, ApJ, 582, 1036 
Nuth, J. A., Ill, Wilkinson, G. M., Johnson, N. M., & Dwyer, M. 2006, ApJ, 644, 
1164 

Palacios, A., Charbonnel, C, Talon, S., & Siess, L. 2006, A&A, 453, 261 
Skelton, R.T., Kavanagh, R.W., & Sargood, D.G. 1987, Phys. Rev. C, 35, 45 
Smith, V. v., & Lambert, D. L. 1990, ApJS, 72, 387 
Sweigart, A. V., & Mengel, J. G. 1979, ApJ, 229, 624 

Thielemann, F.-K., Arnould, M., & Tmran, J.W. 1986, in "Advances in Nuclear 
Astrophysics", eds. E. Vangioni-Flam et al., (Gif-sur-Yvette: Editions 
Frontieres), 525 

Vassiliadis, C, & Wood, PR. 1993, ApJ, 413, 641 

Vogelaar, R. B., Mitchell, L. W., Kavanagh, R. W. et al. 1996, Phys. Rev. C, 53, 
1945 

Wasserburg, G. J., Boothroyd, A. I., & Sackmann, I.-J. 1995, ApJL, 447, 37 
Zinner, E. 1998, Annu. Rev. Earth Planet. Sci., 26, 147 

Zinner, E., Amari, S., Guinness, R. et al. 2007, Geochim. Cosmochim. Acta, 
accepted 

Zinner, E., Nitder, L. R., Alexander, C. M. O.'D., & Gallino, R. 2006a, New 

Astronomy Reviews, 50, 574 
Zinner, E., Nitder, L. R., Galhno, R. et al., 2006b, ApJ, 650, 350 



